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Abstract 

We study the problem of Bose-Einstein condensation in the perfect Bose gas in the 
canonical ensemble, in anisotropically dilated rectangular parallelpipeds (Casimir 
boxes). We prove that in the canonical ensemble for these anisotropic boxes there 
is the same type of generalized Bose-Einstein condensation as in the grand-canonical 
ensemble for the equivalent geometry. However the amount of condensate in the indi- 
vidual states is different in some cases and so are the fluctuations. 
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1 Introduction 



Many calculations in the grand-canonical ensemble (GCE) show a dependence of Bose- 
Einstein condensation (BEC) on the way the infinite volume limit is taken. For example, in 
P and |3] the authors study the the perfect boson gas (PBG) in the GCE in rectangular 
parallelepipeds whose edges go to infinity at different rates (Casimir boxes, see [?f). They 
showed that this anisotropic dilation can modify the standard ground-state BEC, converting 
it into a generalized BEC of type II or III. For a short history of the notion of generalized 
BEC we refer the reader to [T] and jHj- 

On the other hand, due to the lack of (strong) equivalence of ensembles, the PBG in the 
canonical ensemble (CE) and in the GCE gives different expectations and fluctuations for 
many observables. For example, it was shown in [5 4 that for the isotropic dilation of the 
canonical PBG the distribution of ground-state occupation number is different from the 
one in the GCE. The same is true for the fluctuations of the occupation numbers, which 
are shape dependent and are not normal or Gaussian. Therefore this lack of equivalence of 
ensembles does not allow us to deduce the same shape dependence for BEC in the CE as for 
its grand-canonical counterpart and so far the question of whether it is true in the CE has 
not been considered except in a special case p. 

The aim of the present paper is to fill this gap. We study the problem of BEC in the PBG 
in the CE, in anisotropically dilated rectangular boxes. We shall prove that in the CE for 
these anisotropic boxes there is the same type of generalized BEC as in the GCE for the 
equivalent geometry. However the amount of condensate in the individual states is different 
in some cases and so are the fluctuations. 

We would like to note that there is a renewed interest in generalized BEC both from the 
theoretical jEj, [12], [H] and experimental (THj, [E] point of view. This due to recent exper- 
iments which produce "fragmentation" of BEC (see e.g. jH], [E]), that is, the single state 
condensation can be "smeared out" over two or more quantum states. We return to this 
point in Section 4. 

The structure of the present paper is as follows. In the rest of this section we give the 
mathematical setting. In Section 2 we collect together the results about PBG in the GCE 
that we shall need. In Section 3 we study the PBG in the CE for the system of anisotropic 
parallelpipeds. We start by giving some results which are common to the three cases cor- 
responding to the three characteristic ways of taking the thermodynamic limit. These are 
determined by how fast the longest edge grows: (a) faster than the square root of the volume, 
(b) like the square root of the volume and (c) slower than the square root of the volume. In 
the three subsections of Section 3 we study these cases separately. In Section 4 we discuss 
the results. 

We finish this section by establishing the general setting and notation. 
Let A v be a rectangular parallelepiped of volume V : 

A v := {x e M 3 : < x j < V a >,j = 1,2,3} , (1.1) 

where 

«i > «2 > «3 > 0, and «i + a 2 + «3 = 1. (1.2) 



The space of one-particle wave-functions is H Y = L 2 (A V ) and the one-particle Hamiltonian 
t v is the self-adjoint extension of the operator —A/2 determined by the Dirichlet boundary 
conditions on dA v . We denote by {Ek(V)}^ =1 the ordered eigenvalues of t v : 

< E X (V) < E 2 {V) < E 3 (V) < . . . . 

We also introduce the boson Fock space on Ti v defined by J r (H v ) := ©™ =0 ^" symm , where 
:= \(& n j= i'H v \ stands for the space of n-particle symmetric functions. Then 

V / symm 

Ty denotes the n-particle free Hamiltonian determined by t v = Ty 1 ^ on 7i™ symm , and T v 
the corresponding Hamiltonian in the Fock space. 

Now the expectations for the PBG in the canonical ensemble at temperature and density 
p = n/V are defined by the Gibbs state 



(-f v (p) := (ZyinT'Trnn^ (-) C -/*V , (1.3) 

where 



Z v (n):=Tr n , sy _e-< n) (1.4) 

is the n-particle canonical partition function. As usual we put Z v (0) = 1. The grand- 
canonical Gibbs state is defined by 

(-f v C (aO := (Sv^r 1 Tr^ (Wv) (-) e -W*-i*v) t (1.5) 

where /i is the corresponding chemical potential. Here N v is the particle number operator, 
that is, N v := Ylk>i where denotes the operator for the number of particles in the k-th 
one-particle state. The grand-canonical partition function at chemical potential fi < E\{V) 
is 

S v (u):=Tr^ Hv) e-^-^. (1.6) 

Because of their commutative nature it is useful to think of N v and {Nk} k>1 as random 
variables rather than operators. 

Notice that the one-particle Hamiltonian spectrum a(t v ) = {Ek{V)}'^ =1 coincides with the 
set 

( 7T 2 3 n 2 ] 

e„y = Y ^^-:^ = 1,2,3,...;, j = 1,2,3|, (1.7) 

described by the multi-index n = (71%, n 2 , n 3 ). Then and the ground-state eigenvalue Ei(V) = 

e (l,l,l),V- 

Let {r)k(V) := E k (V) - E^V)}^. For a given V we define F v : R -»■ R + by 

Fv(v)--=y#{k:Vk(V)< V }. (1.8) 

Note that Fy is a nondecreasing function on R with F v (rj) = for r] < 0. VF v (i] — Ei(V)) 
is the distribution of the eigenvalues (integrated density of states) of the one-particle Hamil- 
tonian t Y . One can prove in many ways, for example by using Lemma f3. II or by taking the 
Laplace transform, that 

F(rj) := lim F v (r}) = (V2/3n 2 )r] 3/2 , 7] > 0. (1.9) 

V^oo 



We shall show (see Lemma EZQ) that F v (r}) < (V2/3n 2 )rj 3 / 2 if 77 > C/V 2 " 3 for some C > 0. 
This bound and ()1.9|) imply that the critical density of the PBG: 



pc ■= J ^rrr F ^ = / 7^i F{d7]) < °°- (L10) 

(e,oo) 

is finite for any non-zero temperature. 

By (|1.5jl the mean occupation number of the PBG in the grand-canonical ensemble in the 
state k is given by 

Let p v (p) < Ei(V) be the unique root of the equation 

P=^(N v f v c (p) (1.12) 

for a given 1/. Then a standard result |3] shows that the boundedness of the critical density 
(jl.lOj) implies the existence of generalized BEC with condensate density, p , given by: 

p := lim Jim — V (iV fe ) s v c (p v (p)) = P - Pc for p > p c . (1.13) 

eiO V—*oo V ^— ' 

{fc:£ fc (V)<e} 

Following the van den Berg- Lewis- Pule classification £Q and [3] , it is useful to identify three 
categories of generalized BEC: 

I. The condensation is of type I when a finite number of single-particle states are macro- 
scopically occupied. 

II. It is of type II when an infinite number of states are macroscopically occupied. 

III. It is of type III when none of the states is macroscopically occupied. 

For a specific geometry we have more detailed information at our disposal. In the next 
section we collect the results from pQ that we shall need later about the GCE in the case of 
the anisotropically dilated parallelepipeds ([Lip . 

Remark 1.1 Though we have chosen here to work with Dirichlet boundary conditions, the 
proofs in this paper can be adapted without difficulty to periodic or Neumann boundary 
conditions. 

Remark 1.2 Note that according to the classification presented above the condensate 
"fragmentation" is nothing but a generalized BEC of type I or II. 



2 Generalized Bose-Einstein Condensation of the Perfect Bose Gas in the Grand 
Canonical Ensemble 

Proposition 2.1 ( 1 , Theorem 1) 

Let p v (p) = Pv(p) — E\(V). Then the behaviour o//Z v (p) is as follows: 

1. For p < p c , limy^oo Pv(p) = p(p) where p(p) < is the unique root p(p) < of the 
equation 

00 





2. For p > p c , liniy^oo fi v (p) = and for V — > oo ; 

Mp) = - {PV(p - Pc )Y l + 0(1/1/), if an < 1/2; (2.2) 
7V(P) = - {PVA(p)}- 1 + 0{1/V), if ax = 1/2; (2.3) 
TLM = - {2(3V 2 ^\p - p,) 2 }" 1 + 0(1/V), if ai > 1/2, (2.4) 

where A(p) is the unique root of the equation 



(P ~ Pc) = Yl 

3=1 



TV 2 



(f-l)+ A- 1 



-i 



(2.5) 



The next statement by the same authors shows that there are different types of generalized 
BEC corresponding to different asymptotics ()2.2|) - (j2.4j) . 

Proposition 2.2 (pQ) 

For p < p c there is no generalized BEC and therefore no BEC of any type. 
For p > p c there is generalized BEC and all three types of BEC occur: 

1. For ai < 1/2 only the ground-state is macro scopically occupied (BEC of type I): 

lim i {NX (Mp)) = ( r Pc ' f r n Z n !' u ( 2 - 6 ) 

2. For oi\ = 1/2 there is macroscopic occupation of an infinite number of low-lying levels 
(BEC of type II): 

lim 1 W v ° MP)) = { { n {nl ~ 1)1X2,2 + A " Yl ' f n : h !' }]' (2.7) 
y-+oo V /v v /y [0, for n 7^ (ni, 1, 1). 

5. Finally, for ai > 1/2 no single-particle state is macro scopically occupied (BEC of type 
III): 

lim y- 1 (iV n ) 6 v c ( / i v -(p)) = 0, (2.8) 
lim V^-V (N D )f (MP)) = ( ' (P " Pc)2 ' f n : h' Ml' (2-9) 

V^oo IV \r \rj, y /Of II f («!, 1, 1). 

We shall need an easy generalization of the foregoing proposition to obtain the distribution 
of the random variables through their Laplace transform. 

Theorem 2.1 Let p > p c . Then: 

1. Fora x < 1/2, 

hm /expf-A^V(Mp)) = | 1 + Kp-PeV f° rn =^ 1 ^ (2.10) 



V-+00 \ \ V 



v 



0, for (1,1,1). 



2. Fora x = 1/2, 

I f AT \ \ sc { {n 2 1 -l)ir 2 /2 + A- 1 

x v J 1 v { 0, for n ^ (m, 1, 1). 

(2.11) 

5. For ai > 1/2, 

Mm ^exp (-A^ ^ (A*v(p)) = 1, (2.12) 



V 



r / ( \ N " \\f J ^777 V> fom=(n h l,l), 

\ GXP ( ~ X V^ ) ) (Mp)) = 1 + 2A (p - Pc ) 
\ v / / v [0, /orn/ (ni,l,l). 

(2-13) 

Proof: This follows easily from Proposition 12.11 and the identity: 

1 _ e -0(»7k(V)-FvOO) 

(exp (-Atf*)> fl v ° Wrf) = 1 _ e .^ (v0 .p v(p)+V/?) (2-14) 

□ 



We shall require some properties of the Kac distribution K\(p;dp), see e.g. [U EH EU EH] • 
The Kac distribution relates the canonical (jl.3|) and grand-canonical (jl.5j) expectations in a 
finite volume: 

(-f v c (p) = j (-f v (x)K v (v,dx). (2.15) 

[0,oo) 

The limiting Kac distribution gives the decomposition of the limiting grand-canonical state 
(— ) 5C (p) into limiting canonical states (— ) c (p). In the particular case of the PBG it is more 
convenient to define the Kac distribution in terms of the mean particle density, rather then 
the chemical potential. Therefore we define 

K v (p ; dx) := K v (p v (p); dx) (2.16) 

so that 

Rv C (Mp))= J (-) C v (x)K v (p;dx). (2.17) 

[0,oo) 

The next proposition proved in pP gives the limiting Kac density for anisotropically dilated 
parallelepipeds: 

Proposition 2.3 Let 

K(p;dx) : = lim K v (p;dx). (2.18) 

V^oo 

If P < Pa then the PBG limiting Kac distribution has the one-point support: 

K(p;dx) = 5 p (dx). (2.19) 

If p > p c , then: 



1. For a x < 1/2, 



K(p ; dx) 



1 / x - p c 
exp 

P ~ Pc \ P- Pc 



for x < pc, 
dx, for x > p c . 



2. Forct! = 1/2, 



K(p ; dx) = < 



0, 

7T 2 sinh(2/A - tt) 
(2/A-ir)* 



(2.20) 



for x < p c 



x 



v n=l 



J^(-l) n l n 2 exp ^- (x - p c )y ^n 2 + - ljj dx, for x > p c . 



(2.21) 



3. Fora x > 1/2, 



K(p ; dx) = 5 p (dx). 



(2.22) 



3 Generalized Bose-Einstein Condensation and Fluctuations of the Perfect Bose 
Gas in the Canonical Ensemble 

In this section we prove results for the CE analogous to those for the GCE. We are forced to 
use different methods for the three regimes, so we treat them in separate subsections. But 
first we give some results which will be useful in all three cases. The basic identity for the 
canonical expectations at density p = n/V of the occupation numbers is (see equation 
(10)): 

n 

(exp(-AiV fe )) c v (p) = e A - (e A - l)Z v {n)- 1 £ e~^ E ^ n -^ Z v {m). (3.1) 

m=0 

The canonical expectations are notoriously difficult to calculate and are only accessible 
through the grand-canonical expectations. In the two cases a.\ < 1 and a.i — 1 we shall 
exploit the fact that the sum on the righthand side of equation (|3.1|) is very similar to the 
grand-canonical partition function. 

The next theorem shows that the canonical expectations are monotonic increasing in the 
density. Note that this theorem holds for the PBG with any one-particle spectrum. 

Theorem 3.1 For fixed k > 1 and fixed V, the canonical expectations for the PBG, 
(exp(—XNk)) c v (p), are monotonic decreasing functions of the density p for A > while the 
moments (iV£) v (p), r > 1 are monotonic increasing functions of the density. 

Proof: From ()3.1|) we get 

(exp(-XN k )) c v ((n + l)/V) - (exp(-AiV fc )) c v (n/V) 

= -(e x - l)\e-^ +x ^ n+1 ^Z v (n + l)' 1 



f ^ G 1 Z v (n + 1) Z Y (n\ ) I 



m=0 




Since 

'Z v (m + 1) Z v {m)\ Z v (m) fZ v (m + l) Z v (n + 1) 



Z v (n + 1) Z v (n) J Z v (n + 1)\ Z v (m) Z V (N) 
by the inequalities (see 



Z v (m + 1) Z v (m + 2) Z v (n + 1) 



Z v (m) Z v (m + 1) ^v(") 

and by (|3.2J) we get the monotonicity: 

(exp(-\N k )) c v ((n + l)/V) - (exp(-AiV fc )) c v (n/V) < 0. (3.3) 

By differentiating ()3.2j) r times with respect to A at A = 0, 

W((n + 1)/V)-W(n/V) 

= {n r - (n - l) r }e- /3jB ^ n+1 )Z v (n + l)" 1 

£- n U j ( Jl Uv(n+1) Zy{n))- 



(3.4) 

\ /, , i i i /,>(//) / 

m=0 

which is positive by the same argument. □ 

Remark: In this paper whenever we take the limit 

lim (-)» (3.5) 

V — i-oo 

we shall mean that we take the system with n particles in a container of volume V n = np 
and then let n — > oo, that is, 

Hm (-f np (p). (3.6) 
The next theorem is valid for containers of any geometry and not just for rectangular boxes. 



Theorem 3.2 For p > p c the generalized condensate in the CE at density p is equal to 
p — p c , that is 

hm lim (Nk/Vf v (p)=p- Pc . (3.7) 

Vk<e 

Proof: The statement is true for the imperfect (mean- field) Bose gas in the GCE, see [2]. 
Since the mean-field term in the CE is irrelevant, the theorem follows from monotonicity 
and the fact that the Kac density for the imperfect Bose gas has one-point support. □ 



In the next theorem we shall make certain assumptions that are clearly satisfied for the 
parallelpipeds we are considering. We believe that in fact they hold much more generally. 



Theorem 3.3 Suppose that limy_ >0 o (Nk/V) 5 y(p:A{p')) and K(p'; [p, oo)) are continuous in 
p' at p and that K(p; [p, oo)) ^ 0. Then limy-^ (Nk/V) 5 y (/iy(p)) = implies that 
hm v ^ oo (N k /V) c v (p) = 0. 



Proof: Using the decomposition (J2.17)) and monotonicity we get for any e > 0: 



lim (N k /V) s v c (p A (p + e)) = lim / ( N k /V) c v (x) K v (p + e; dx) 

V^oo V— >oo j 

[0,00) 



> lim / (N k /Vf v (x)K v (p + e;dx) 

V^oo J 
[p,oo) 



> lim sup {N k /Vf v (p) K v (p + e; [p, 00)) 



lim sup (N k /Vf v {p) K(p + e; [p, 00)). 



Since limy^oo (N k /V) v (p A (p')) and K(p'; [p, 00)) are continuous in p', letting £ tend to zero, 
we get 

lim (N k /V) g v c (p A (p)) > lim sup (N k /V)\{p) K(p; [p, 00)) 
v^oo y^oo 

and because K(p; [p, 00)) does not vanish the result follows. □ 



Remark: Note that this lemma implies that for p < p c , there is never BEC in the CE. 

Before looking at the three cases ai < 1/2, ai = 1/2 and a% > 1/2 we first obtain lower and 
upper bounds on the density of states. 



^ - Cy--) 3 < Fv fo) < ^ + El (V)f 2 (3.8) 



Lemma 3.1 

for some C . 
Proof: 

VFyfa - E 1 (V)) = #| n |nGN 3 , y ^^-<r / |, (3.9) 
that is VF v {j] — Ei(V)) is the number of points of N 3 inside the ellipsoid 

%2 1 y2 1 z<1 - 1 fQ mi 

2F 2a ir//7r 2 2V* a *r)/7r 2 2V 2 ^r l /n 2 ' 1 ' ; 

If we associate the point n with the volume of the unit cube centered at n — (|, §, 5) we 
see that this number is less the volume of the ellipsoid in the first octant which is equal to 

n(2V 2ai r]/TT 2 )^(2V 2a2 7]/n 2 )^(2V 2a3 7]/TT 2 )^ /Q = V2t] 3/2 V/3tt 2 . Thus 

VF v (r,-E 1 (V))<^r, 3 / 2 V (3.11) 

and so 

F v ( V )<^(v + E 1 (V)f 2 . (3.12) 



Let a > b > c > 0, let A = 1 — 3/c and a' = Xa, b' = Xb, and d = Ac. If the point in the 
first quadrant (x,y,z) satisfies x 2 /a' 2 + y 2 /b' 2 + z 2 /d 2 < 1, then it satisfies (x + I) 2 /a 2 + 
(y + l) 2 /b 2 + (z + l) 2 /c 2 < 1. That is, each point inside the first quadrant of the ellipsoid 
x 2 /a' 2 + y 2 /b' 2 + z 2 jd 2 = 1 lies in a unit cube with the corner n G N 3 (with n\ > x, ri2 > y 
and n 3 > y) inside the ellipsoid x 2 /a 2 + y 2 /b 2 + z 2 /c 2 = 1. Therefore 

VF v{ ,-E l( V))>^ V (^--^) (3,3) 

yielding 

Fv{v) > ^ (,, + El{Vyr _ > gy ^ _ _|_) 3 . (3,4) 

□ 



3.1 Case ai > 1/2. 

We study this case first because it is the simplest since the limiting Kac distribution is a delta 
measure concentrated at p and we have strong equivalence of ensembles (see Proposition ^. 3j) . 
We shall use this fact together with the monotonicity properties of Theorem 13 .11 to show that 
in this case the limiting canonical and grand-canonical expectations are identical. 



Lemma 3.2 For a,\ > 1/2 and A > the following inequalities hold 

^ vJ^o ( exp (- X V^)))l Mp ~ £)) 
>limsup^exp (-A^^^(p) 

> liminf ^exp (-A^)^y))(p) 

> limsup ton ^exp ^-A ^^y j J (fi A (p + e)). (3.15) 



Proof: We start with the first inequality. Using the decomposition (J2.17)) we get for any 

e > 0: 



hm /exp f-\^-\ \ (p v (p - e)) 

[0,00) 

[0,P) 

> lim sup j ^exp ^-A ^^y j ^ (p) K v (p - e; [0, p) 
= limsup^exp (-A^^(p). 



(3.16) 



In the penultimate inequality equality we have used the monotonicity established in Theorem 
13.11 and in the last one we have used 1)2.19)1 and ()2.22jl . The last inequality in 1)3.15)1 is proved 
similarly: 



= j ^exp f-A ^j^y ^ \ (x) K v (p + e; dc) 
[0,00) 

= y (| ex P f-A^i~^))^ (x) K v (p + e; dx) 

[0,P) 

+ y (exp r-A ^^y j \ (x) K v (p + e; dx) 
[p, 00) 

< K v (p + e; [0, p)) + ^exp ^-X^L-^j ^ (p)K v (p + e; [p, 00)). 



(3.17) 



Therefore 



Km ^exp ("A^t^) ) JPv(p + e)) 

<l^f^«p(-A^^)) C (p)- (3-18) 



□ 



The following theorem and corollary give the distribution and the mean of N n /V 2 ^~ ai>) and 
therefore they give the fluctuations about the mean. 



Theorem 3.4 If a\ > 1/2 and p > p c then the limiting distribution in the canonical ensem- 
ble of iV n /l/ 2 ( 1_Ql ) has Laplace transform for A > given by 

hm (expf-A^^l )(p)= 1 + 2A (p - p c ) 

x v 7 ' ^ 0, for (m, 1, 1). 



(3.19) 

Proof: This follows from the preceding lemma and Theorem 12.11 □ 
Corollary 3.1 If oi\ > 1/2 and p > p c then 



Hm / ^ \' f 2A (p-p c ) 2 , /orn=KM), 

\ ^(l-m) / W \ o, /or n ^ (m, 1, 1). 



Proof: Again we have to check that there exists K < oo such that for all V 



(3.20) 



N, 



k 



y2(l-ai) 



2 



c 



(p) < K (3.21) 



Then the corollary follows from the preceding theorem. We have again for any s > 0: 

1 v (p v (p + e)) 



y2(l-ai) 



v 

2^ 



k 



[0,oo) 4 7 s 

2 



K v (p + e; dx) 



7 / { ( tt£^) ) (*) &v(p + dx) 



y2(l-ai) 



[p,oo) 



V 



>( (t^T) ) ) (p)Kv(p + e;[p,oo)) 



(3.22) 



if is large enough. This implies the existence of K as above since 

(p v (p + e)) 



iV fc 



y2{i- ai ) 



V 



(3.23) 



converges as V — > oo. 



□ 



Corollary 3.2 When ci\ > 1/2, there is type III BEC in the canonical ensemble. 



Proof: From the preceding corollary or from Theorem 13.31 we can deduce immediately that 
for any p > p c . □ 



3.2 Case a x = 1/2. 

For this case BEC into the ground state is treated in Here we extend the result to higher 
levels. 

Because for ai = 1/2 the spectral series (jl.7j) corresponding to n\ = 1, 2, 3, ... has the smallest 
energy spacing tt 2 /2V, it plays a specific role in calculations of the limiting occupation 
densities. Let 

e n := 7r 2 n 2 /2 , rj m>n := (3{e m - e n ) , b m , n := v^ n J| (1 - r] m ,„/?7 m ' in ) _i (3.25) 
for m ^ n. 

In ([S|) the following result was proved: 
Let cti = 1/2. Then for p > p c 

oo 

^ b m ,i {Vm,i(p ~ Pc) ~ 1 + exp [-rj m>1 (p - p c ))} 

oo 

^2 b m,iVm,i {1 - exp [-r) m>1 (p - p c ))} 



771=2 



Here we give an extension of ()3.2(ij) to other fc's. Note that by Theorem 13 .31 and comparison 
with the GCE, in this case there can only be condensation in states corresponding to n = 
(ni, 1, 1). The main tool in the technique developed in [S] and jB] is the following identity: 

Let {Kk )V {dx)} k>1 be (non-normalized) measures whose distributions are the functions 

K M-S Z ^ r ) ex P i-P(VPk ~ rE k (V)} , for r/V < x < (r + 1)/V, ( . 
%W "{0, forx<0, {6 - Z(} 

for r = 0, 1, 2, . . . and some {pk\k>\- Then we can re-write equation (J3.1)) as follows 

(exp{-\N k /V}} c v (p) = e~ Xp je Xx K Kv (dx)/K Kv (p+l/V) 

[O.P+l/V] 

rp+l/v 

= e- x / v -\e~ xp K kjV (x)e Xx dx/K kjV (p+l/V). (3.28) 
Jo 



We shall use this identity to calculate the thermodynamic limit of its left-hand side for a 
given density and k > 1. From ()3.27|) we can calculate the Laplace transformation of the 
measure K k)V (dx): 

J e~ Xx K Kv (dx) = (1 - e - x ' v )e- v ^ E v (E k (V) - X//3V) (3.29) 

k 

where for the PBG the grand-canonical partition function (jl.fi |) has the explicit form: 

oo 
k=l 

Now we fix the p k s which are still arbitrary by the defining 

Pk .- _ J_ J^i n |i _ e -m<y)-*.0O)\ (3.30) 

and prove the following lemma. Let 

K n ,v '■= ^(n,l,l),V (3.31) 



Lemma 3.3 Let aj = 1/2. Then 
K n (x) := lim K n v (x) 



(3.32) 



/or x < p c 



.!)("!) ^ 6m,»»7tn,»{l-exp[-»7n*,n(^-/3c)]} for X > p c . 
m=l, mj^n 



Proof: Let A > V(E k — J5i). From the definitions flQSD and (|3~5Uj) we get: 



e- Xx K kjV (dx) = Y[- 



1 _ g-jSC^oo-SfcOO)! 



exp < — In 



e -/3(i? J (V)-i? fc (V)+A//3y) 
j^ft k j+k 



1 + 



1 - e- A /v 



;/ 3(£j(V)-E fc (K)) _ J 



For > 1 and 77 > 771 (V) — rj k (V) we define the shifted integrated density of states, cf (jl.8j) 
FkAv > 0) := ^# {J : %W < V + Vk(V), J ? k} = F V ( V + Vk (V)) - ^l [0 ,oo)(v), (3-33) 



For b > a, let 



I kjV (a,b) :=V J In 

[a,6) 



1 - e-W 



el* - 1 



F k ,v{dr]) 



(3.34) 



Then 



2> 



1 - e- x ' v 



J3{E 3 (V)-E k (V)) _ l 



= V j In 

foi(V)-»jkCV),oo) 



1 + 



1 - e- A / y 



e /*7 - 1 

i k , v (m(v)-Vk(v),oo). 



F k , v {drj) 



(3.35) 



We can write 

Ik,v(Vi(V)- V k(V),oo) = h,v(Vi(V) -Vk(V),l/V 2a3 ) +I k , v (l/V 2a3 ,oo). (3.36) 
Since limy^oo rj k (V) = 0, by Lemma I3~T1 we have 



Jim F k;V {rj) = Km F V ( V + n k (V)) = F( V ) 

V— >oo K^oo 



and 



f k M < c'v 3/2 

for rj > 1/V 2a3 , using the estimate x — x 2 /2 < ln(l + x) < x we get 

r2a 3 



Let 
Then 



lim I k>v (l/V' a3 ,oo) = \p c 



G k>v (0 := VF k J£/V). 



hAvi(V)-Vk(V)A/V 2a3 )= j In 



1 + 



1 - e- x ' v 



e K/v _ i 



(3.37) 
(3.38) 

(3.39) 
(3.40) 

G k , v (d£). (3.41) 



Now let E k (V) correspond to £( n ,i,i),v Then in the limit, G k ^ v gives a non-trivial point 
measure concentrated on the set {fi~ x r\ m ^ m ^ n}: 



lim G kjV {£ > 0) = # {m : /T 1 



V^oo 



(3.42) 



Therefore, by ((HUD and (jSjlSj) we get 



lim {/ fc)V (^(»7i(V) -%(¥)), 1/F 2c * 3 )} = £ln 



1 + 



<]m,n 



(3.43) 



Thus for A > |?7i ir , 



oo 

lim / K k v (dx)e~ Xx = exp{— Ap c }exp < — / In 
'-><x>J \ J ~ 



V 

o ^ m A 

and the lemma follows by inverting the Laplace transform. 



T)m,n 



(3.44) 



□ 



Theorem 3.5 Let p > p c and a± = 1/2. Let E k (V) correspond to ei n ,i,i),v Then 



V b 



Vm.n 



Km (exp{-AiV fc /^}) c v (p) 



771=1, m^n 



A 



{exp [X(p - p c )\ - exp [-7] m ,n{P - Pc)}} 



^ &m,n?7m,n {1 - exp [-r] m<n (p - p c )}} 
m=l, m^n 



(3.45) 



Proof: The identity (J3.28)) gives: 

p 

Xe' xp [ K k (x)e Xx dx 



lim (exp {-\N k /V}f v (p) = 1 °— . (3.46) 



From the preceding lemma and for A > 

p 



Xe~ xp I K k (x)e Xx dx & m,„W Ae ~ V J e Xx {1 - exp [-7? m , n (x - p c )\} dx 

= ^ Pc 

K k(p) ^ b m, n r]m,n {1 ~ exp [-p m ,n(p - Pc)}} 



^ &m,n {p m ,„exp [-A(p - p c )] - A exp [-r/ m ,„(p - p c )] - (r) m , n - A)} 

, Vm,n 

b m,nVm,n {1 ~ exp [-7? m ,„(p - p c )]} 

(3.47) 

This gives (l3~43|) . □ 



We are now in a position to prove that in this case there is BEC of type II for p > p c . 

Theorem 3.6 For p > p c and for a% = 1/2 all states with n = (n, 1, 1) are macro scopically 
occupied (BEC of type II) while all the other states are not. The occupation density is given 
by 



^ b min {p m ,„(p - p c ) - 1 + exp [-r) m ,n(p - Pc)]} 

,. / AW,1) \ / \ m=l,^n 

^{v/J p) = 

2 , &m,n^m,n {1 - exp [-?7 m>n (p - p c )]} 
m=l, mj^n 

Proof: As we mentioned above it is sufficient to check that there exists K < oo such that 
for all 1/ 

<(AV^) 2 »)<^- (3-48) 
Then the theorem follows from the preceding one. We have: 

((N k /V) 2 f v C (p v (p)) = J ((N k /V) 2 ) C v (x)K v (p;dx) 

[0,oo) 

\2\C 



> J ((N k /V) 2 ) v (x)K v (p;dx) 

[p,oo) 

> ((N k /V) 2 ) C v (p)K v (p;[p,oo)). 

(3.49) 

This implies the existence of K as above since {{N k /V) 2 'f^ (p v (p)) converges as V — > oo 
and K v (p; [p, oo)) converges to a non-zero limit. □ 



3.3 Case ai< 1/2. 



In jHj the canonical PBG in parallelepipeds 

A v := {x£K 3 :0< x s < ajV 1/3 ,j = 1,2,3} , ai a 2 a 3 = 1 (3.50) 

were considered. It was proved that for this system there is BEC of type I . In particular it 
was proved that: 

Proposition 3.1 ( 5 , Theorem 1) For the PBG in parallelepipeds A3.5U\) . the following limits 
hold when A G R: 

if p > Pc, and 

lim (exp {-\vm k /V}) c v (p) = 1 (3.52) 

V— >oo 

for any < 7 < 1 and k > 1, if p < p c . 

Note that in Theorem 1, A < but this is not necessary. 

In this section we shall prove a similar result to Proposition 13.11 for case of the rectangular 
parallelepipeds (jl.lj) with a.\ < 1/2, that is, we shall show that in this case there is also 
type I BEC. It is sufficient to show that there is condensation in the ground state since by 
Theorem 13.31 no other state can be macroscopically occupied. 

Let Ki <v (dx) be as in ()3.27|) in Section 3.2. 
Lemma 3.4 Let a\ < 1/2. Then 

lim K ljV (dx) = 5 Pc (dx). (3.53) 

Proof: The proof is almost identical to that of Lemma 13.31 The only difference is that since 
Vj(y) ~ Vi(Y) > aj,k/V 2a \ (j3~?nj) implies that 

G hv (C > 0) = # : r)j(V) - Vl (V) < i/V} -> when V - 00, (3.54) 

for oil < 1/2 and the lemma follows. □ 

We can now prove that in this case there is BEC of type I for p > p c . 

Theorem 3.7 For p > p c and for ai < 1/2 only the ground-state is macroscopically occupied 
(BEC of type I): 

v->co V [ 0, /or nyt (1,1,1). v ' 



Proof: From the preceding lemma and the identity ()3.28j) we have for A > 

Jim (exp (-XN a /V)f v c (p) = -X(p - p c ). (3.56) 

It is sufficient to show the second moment is bounded, that is, there exists K < oo such that 
for all V 

((N k /V) 2 f v (p) < K. (3.57) 
The bound can be obtained by the same argument as in Theorem 13.61 □ 



In the rest of this subsection we study the fluctuations of N\/V. We need the shifted 
integrated density of states in d dimensions d = 1,2, 3, in the unit boxes [0, l} d : 

F} d \v) ■= # |n | n G N d , y $>i " l ? < ^} > ( 3 - 58 ) 



Theorem 3.8 Suppose a% < 1/2 and let 7 = 1 — 2a\ > 0. Then for p > p c , 

Jim (exp {XV^N./V - {N x /Vf y {p))}) c v {p) = 

exp (gi(\)) , if a 3 < a 2 < oix < 1/2, 

exp (2</i(A) + fla(A)) , if a 3 < a 2 = «i < 1/2, 

exp (3pi(A) + 3g 2 (X) + g 3 (X)) , if a 3 = a 2 = a x = 1/3, 

(3.59) 

where 



9iW 



»(A) 



<? 3 (A) 



(0,00) 



ln ' 1 + Pv) + Pv 



F^(drj), 



ln 1 



A 



A 



(0,oo) 2 



In 1 + 



P(Vl+V2)J P(Vl + V2) 

A 



+ 



F^\d Vl ,d V2 ) 
A 



(0,oo) 3 



P(r]i +V2 + V3)/ P(Vi +^2 + 773). 



Fr\d Vl ,drj 2 ,drj 3 ). 



(3.60) 



Remark: Note that 3pi(A) + 3g 2 (A) + (73(A) is the same as g(X) in ( e 5j). 
Proof: Let 

L v (x) : = 

Zy(r) exp {-^(Vpi - rE k (V)} , for ^(r/V - p v c )) <x< V^{(r + 1)/V - p, 
0, for x < -V^p v c , 



for r = 0, 1, 2, . . ., where p\ is as in ()3.30|) . 



1 



' e Pv-l 

(0,00) 



F v {dr,) 



(3.61) 



and where we put Z v (0) = 1. Then 

(expiXV^NjV - (p - P :))}f v (p) = J e- x *L v (dx)/K iy (p + l/V) 

(-co, ay) 

where a v = V y (p- p v c ) + V~ 2ai . By Lemma ITU for p > p c , limy^ K iy (p + l/V) 
therefore 



(exp {XV^Ni/V - (p - pl))}f v {p) = Jim e^L v (dx). 

(—00, 00) 



Now 



In J e- Ax L v (dx) = V 
(—00, 00) (0, 00) 



n^(l,l,l) 



In 1 



-In 1 + 



1-e 



1-e 



e Pv - 1 



+ 



+ 



eft - 1 

xv- 2ai 



F v {d V ) 



e P( e n,V— £(1,1,1), v) \ j g/3( e n,V — e(l,l,l),v) \ 

Consider first the case a.3 < «2 < &i < 1/2. We write 

In Je~ Xx L v {dx) = A v +B, 

(—00, 00) 



where 



In 1 + 



1-e 



-\v- 2a i 



+ 



XV 



-2oi 



3^( e (n 1 ,l,l),V- £ (l,l,l),v) — 1 / e 0( € (n 1 ,l,l),V- e {l,l,l),v) _ \ 



and 



(n 2 ,n 3 )^(l,l) 



In 1 + 



1-e 



-\v- 2a i 



+ 



-$(txi,V — v) — 1 / g^( e n,F-£(l,l,l),T/) _ ]_ 



Note that by definition 



A, 



(0,00) 



-In 1 + 



1-e 



+ 



Ay -2ai 



F^(dr,). 



Using the bounds (obtained from the inequality — x < — ln(l + x) < \x — x) 



< 



A - 1 + e 

e x - 1 



< - In 1 



1-e 



-A 



A 



< - 



1 / l-e- A A 2 A-l + e~ A cA 2 



2 V e x - 1 



e x - 1 



e x — 1 j e x — 1 

< 



we get using the Dominated Convergence Theorem 



lim A v = g^X). 

V^oo 



(3.70) 



Using the same inequality 

\2 

< B v < 



E 



4A J 



E 



lT 4 j3 2 V 4ai ({n\-l) , (n|-l) . 

(n2,ri3)^(l,l) I y 2ai 1 \F^2 ' V ia 3 



^ (.,^(1,1) (K " 1) + (*2 - 1)^(— ») + (n| - l)^ (ai -a 2 ))2- 



(3.71) 



Now the summand in the last sum tends to zero as V — > 00 and it is bounded above by 

1 



(K - 1) + (n| - 1) + (n§ - 1)) 



2 • 



Since 



E 



1 



i „,„:„,m.i, (K-l) + (nI-l) + (n|-l)r 
B v — > as V — > 00 by the same theorem. 
Next we consider the case 03 < 0:2 = cm < 1/2. Now we take 



< 00, 



(3.72) 



(3.73) 



A, 



and 



E 



(ni,n 2 )^(l,l) 



-In 1 + 



1-e 



-Ay- 2 «i 



+ 



3 /^( e (?l 1 ,?l2,l),V~ e (l,l,l),V r ) ]/ g/5( <E (n 1 ,n 2 > 1 ). v ' — e (l,l,l),v) \ 



(3.74) 



71 3 ^l 



In 1 



1-e 



3^(e n ,V— e (l,i,l), v) — 1 / g/3(en,V— e(i,i,i),v) _ ]_ 



(3.75) 



In this case by definition ()3.58|) 



.4, 



2£ 

ni^l 



In 1 + 



1-e 



+ 



J3{^(nx,\,l),V~ € i\,l,l),v) \ J Q0( e (n 1 ,l,l),V~ e (l,l,l),vO J 



+ E 

rai^l,n 2 ^l 



In 1 + 



1 - e 



3 /3( e (n 1 ,n 2 ,l),V — e (l,l,l),v) I g/3( e (n 1 ,n 2 ,l),V r— e (l,l.l),v) \ 



(0,00) 



/ 1 _ e -AV- 2 "i \ X y-2 ai 
+ e /3y- 2o! i77 _ ]_ J + e /3V- 2a i V _ 1 



+ 



(0,00) 



In 1 + 



1-e" 



+ 



XV 



-2ai 



oj3V- 2a i{ Vl + m ) - 1 ePV-^iim+m) _ 1 



Fl 2 \drn,dr]2). 



(3.76) 



By the same argument as above 



lim A v = 20i (A) + 02(A) 

V— >oo 



and 



lim B v = 0. 

V->oo 



(3.77) 



(3.78) 



Finally for a 3 = a 2 = ot\ = 1/3 



L v (dx)e 



-Xx 



(—00, 00) 



3£ 

m^i L 



In 1 + 



1-e" 



+ 



+3 £ 



-In 1 + 



e /3( e (m 1 ,l,l),V' — e (l,l,l),vO — 1 / e^( e (™l,l,l),^ _(E (l,l,l),v) _ ]_ 

AK- 2q i \ ^y-2ai 



1 - e" 



3 /3(e(n 1 ,n 2 ,l),V — e (l,l,l),Vj ]_ I g^( e (n 1 ,n 2 ,l),V — e (l,l,l),v) ]_ 



(0,00) 



+ 3 



(0,oo)2 



+ 



(0, oo) 3 



In 1 + 



In 1 + 



In 1 + 



1 - e ~ xv ~ 2ai \ XV~ 2 ^ 



+ 



e [3V- 2a iri — \j e /3V- 2a lr) _ I 



1-e 



XV~ 2a ' 



e /3V- 2a i(m+V2) — I ) e pv- 2a i( ril + m ) _ 1 



F?\d Vl ,d V2 ) 



1-e 



\y-2ai 



e pv- 2oi i(rn+V2+ri3) - 1 / ePV-^iirn+m+vs) _ 1 



(3.79) 



Thus 



lim In \e- Xx L v (dx) = Zg x {\) + 30 2 (A) + g 3 (X). 



(3.80) 



(—00, 00) 

We have therefore proved that 

lim (exp {XV^NJV - (p - p^)}% (p) 

V — >oo 



exp (01(A)) , if a 3 < a 2 < a x < 1/2, 

exp (2#i (A) + 02(A)) , if a 3 < a 2 = «i < 1/2, 

exp (30i(A) + 30 2 (A) + 03(A)) , if a 3 = a 2 = aci = 1/3. 



To finish the proof we centre the distribution about (Ni/V) c v (p). From ()3.62|) we get 
({XV^NJV -(p- pl))f) C v (p) = J x 2 L v (dx)/K iy (p + 1/V) 

(—00, ay) 

< 2 Jx 2 L v {dx). 



(—00, 00) 



Using (|H.fj4jl this gives 

(—00,00) 
3 yl-4ai f l 



< ^V(^). (3.81) 

(—00,00) 

Thus by the same arguments as above this second moment is bounded. Since g[(0) = g' 2 {0) = 
#3(0) we then have 

lim V^m/Vf, (p) - (p - p%) = lim (V 7 (Ni/V - (p - p^)% (p) = (3.82) 

V— »oc v— »oo 

completing the proof. □ 



4 Conclusion 

(a) Since the paper by Buffet and Pule [S] it has been known that there are differences in 
the fluctuations of the PBG condensate in the canonical and grand-canonical ensembles. We 
have shown that the picture becomes even more complicated if one passes to the case of 
Casimir boxes where there is already generalized BEC in the GCE. 

We have shown in general that there is a kind of stability principle relating the two ensembles: 
condensation in the GCE is more stable than in the CE. In fact we proved (Theorem 3.3) 
that that the absence of the macroscopically occupied single-particle states in GCE implies 
the same in the CE, whereas the converse is not necessarily true. However in the case of 
the Casimir boxes considered here the two ensembles exhibit the same types of BEC for the 
same geometry. What varies in some cases is the fluctuations and the amount of condensate 
in the individual levels. 

(b) As we have mentioned in Section 1, BEC of type I and II is also known in the literature 
describing the experiments with trapped bosons as "fragmentation" of the condensate, [Oj, 
[TT)j , [B] . In these papers the authors relate this phenomenon to the interaction properties of 
the Bose gas, arguing that it is the exchange interaction that causes bosons with repulsive 
interaction to condense into a single one-particle state whereas for attractive interactions the 
condensate may be "fragmented" into a number of degenerate or nearly degenerate single 
particle states, see Here we have shown that BEC of type II occurs in the non- 
interacting Bose gas and that this is due simply to a geometric anisotropy of the boxes 
known since Casimir On the other hand, there are exactly soluble models (in cubic 
boxes) showing that some truncated repulsive interactions are able to convert BEC in the 
ground state into a generalized condensation of type III (see [TI], [1J). In ^3] an even simpler 
repulsive interaction than in [TT] and jlj is proposed that produces type I condensation in a 
few degenerate single-particle states. 
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